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Abstract. In this talk I discuss the generation of a stochastic background of gravitational 
waves during a first order phase transition. I present simple general arguments which explain 
the main features of the gravitational wave spectrum like the k 3 power law growth on large scales 
and a estimate for the peak amplitude. In the second part I concentrate on the electroweak 
phase transition and argue that the nucleosynthesis bound on its gravitational wave background 
seriously limits seed magnetic fields which may have been generated during this transition. 



1. Introduction 

The Universe has expanded and cooled down from a very hot initial state to (presently) 2.7°K. 
It seems likely that it underwent several phase transitions during its evolution of adiabatic 
expansion. It has already been proposed in the 80ties, that a cosmological phase transition 
can lead to the generation of a stochastic gravitational wave background [1]. If the phase 
transitions are of second order or only a crossover as it has been obtained from standard model 
calculations, we do not expect them to lead to an observable cosmological signal. However, if a 
phase transition is first order, it proceeds via bubble nucleation which is a highly inhomogeneous 
process. Some of the kinetic energy of the rapidly expanding bubbles is transferred to the cosmic 
plasma and leads to turbulence. If small seed magnetic fields are present, they are amplified by 
the turbulent motion which becomes MHD turbulence. Since the fluid and magnetic Reynolds 
numbers in the hot cosmic plasma are extremely high, MHD turbulence persists even after the 
phase transition, where it freely decays. These inhomogeneities int the energy and momentum 
of the cosmic plasma generate a stochastic background of gravitational waves. Here we want to 
investigate this background. 



1.1. Events in the early Universe 

The following important events have most probably taken place in the early Universe and they 
may, under certain circumstances have induced a background of stochastic gravitational waves. 

• Inflation is supposed to have generated the scalar fluctuations in the geometry and energy 
density which have led to the observed anisotropies in the cosmic microwave background 
(CMB) [2, 3]. Simple inflationary models do not only induce scalar perturbations but the 
quantum fluctuations in the gravitational spin-2 degrees of freedom also lead to a scale 
invariant background of gravitational waves of similar amplitude. It is one of the most 



important goals of present CMB experiments like Planck [4] to observe the B-polarization 
which this gravitational wave background should induce in the CMB photons. 
P re-heating: At the end of inflation the inflaton decays into other degrees of freedom 
which eventually reheat the Universe and lead to a thermal bath of relativistic degrees 
of freedom containing especially also the standard model particles. The details of this 
process are complicated and in some cases they generate inhomogeneities leading to a 
significant stochastic background of gravitational waves [5]. Using the standard relation 
between temperature and time in a Friedmann Universe filled with relativistic particles [3] , 
we obtain for a reheat temperature Tj 

Ti ~ 10 14 GeV, t J = 2.3sec(^^) 2 5cff (T)- 1 / 2 ^10- 32 sec, (1) 

* = u \uu\ =2*i(l + *) - 10~ 7 sec, Wi>-~10 7 Hz, (2) 
a(ti)H(ti) f]i 

9cS (T) denotes the effective number of relativistic degrees of freedom at temperature T [3] . 
To find the conformal time rji we have used that in a radiation dominated universe 
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where % = aH is the comoving Hubble rate and a is the scale factor normalized to one 
today, a (to) = 1. Furthermore, since g e ff(aT) 3 is constant during adiabatic expansion 
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so that 

71 = °- 67 x 1010 (iw) " L5 x 10l ° sec (^) 5eff(T) " 1/6 • (3) 

We have set ks = 1 and measure temperature in MeV. In the following we use units with 
c = h = 1 so that length and time have the same units as inverse energy. 
For the >, sign in eq. (2) we have used that the typical frequency generated by such an 
event must be larger than the inverse of the correlation scale which is always smaller than 
the horizon size at the given epoch. Note that Wj denotes the comoving frequency and we 
have normalized the scale factor to unity today, so that co-moving frequencies, length scales 
and times correspond to physical scales today. 

So far no experiment which can detect gravitational waves with frequencies as high as oji 
has been proposed. 

The electroweak transition: In the standard model, for realistic values of the Higgs mass, 
the electroweak transition is not even second order, but only a crossover [6]. However, small 
deviations from the standard model e.g. in the Higgs sector or supersymmetric models [7] 
can can lead to a first order electroweak phase transition at temperature 

T ew ~ 10 2 GeV, t ew ~ l(T 10 sec, (4) 

i] ew ~ 10 5 sec, uj £W > l(T 5 Hz. (5) 

If the electroweak phase transition is strongly first order, the bubbles expand very rapidly 
and the duration of the phase transition is a small fraction, maybe 1% of the Hubble time. 
In this case one expects a peak frequency of w ew ~ W0/r] ew ~ 10~ 3 Hz. This is the frequency 
of best sensitivity for the gravitational wave satellite LISA proposed for launch in 2018 [8]. 



• Confinement transition: Also in this case, standard lattice QCD calculations [9] predict 
a simple crossover. However, the confinement transition can be first order if the 
neutrino chemical potentials are sufficiently large (still within the limits allowed by 
nucleosynthesis) [10]. Models which lead to such high neutrino chemical potentials have 
been proposed for leptogenesis with sterile neutrini [11]. The temperature of the QCD 
transition is about 

T c ~ 10 2 MeV, t c ~ 10~ 5 sec, (6) 
n c ~ 10 7 sec, w c >10~ 7 Hz. (7) 

Gravitational waves in this frequency range could be observed by the pulsar timing array 
which has been proposed recently [12]. 

• Surprises: Last but not least we may hope to find a gravitational wave background from 
a phase transition at some other not predicted temperature which then would indicate 
new physics at the corresponding energy scale, see, e.g [13]. For example, a first order 
phase transition at a critical temperature T* ~ 10 7 GeV may lead to a gravitational wave 
background with a peak frequency of about 100Hz which is the best sensitivity of LIGO. 
As we shall argue below, in a very optimistic scenario it is even conceivable that such a 
background might be detected by the advanced configuration of the LIGO [14] experiment. 

A more comprehensive overview on first order phase transitions as sources for gravitational 
waves can be found in Ref. [15]. 

In the following we present a semi- analytical determination of the gravitational wave (GW) 
signal from a first order phase transition in terms of a few free parameters. We then apply our 
general results to the electroweak phase transition and discuss some consequences. 



2. General results 

Before we outline their derivation, let us present the main results. We consider a phase transition 
at T = T* which generates (relativistic) anisotropic stresses which source gravitational waves. 
We assume that the energy density of this source is px- Since the source is relativistic its 
anisotropic stresses are of the same order of magnitude. The GW spectrum then has the following 
main properties: 

The peak frequency (the inverse of the correlation scale) is larger than the Hubble rate, 
k* ^jH*- Here we always compare comoving quantities which become physical scales today. 

The GW energy density is of the order of 



ttg W (r] ) ~ e S7 rad (??o) f where 

f (H*Arj*) 2 if n,Ar], < 1 , . 

\ 1 if n,Ar], > 1 . 1 J 

Here An* is the duration over which the source is active, typically the duration of the phase 
transition. Q is the density parameter, VLx = px/pc where px is the energy density of the 
component X and p c is the critical density. 

On large scales, k <C k* the spectrum is blue, 



d£lgw{k) ,3 r\ f dk d£lg W (Jz) 

dlog(k) K * ' U9W ~ J T dlog(k) 



The second equality defines the the gravitational wave energy per logarithmic frequency 

interval dU 9w(k) 
interval, dlog(fe) . 



• On small scales, k~^> k* the spectrum decays. The decay law depends on the details of the 
source. If the anisotropic stress spectrum of the source decays with power /c~ 7 , and the 
source is totally coherent in time (see below), the gravitational wave energy density decays 
like k-^ +1 \ If the source is coherent for about one wavelength, the gravitational wave 
energy density decays like /c^ 7-1 ) and if the source is incoherent, the gravitational wave 
energy density decays like k~^^ 3 \ see Ref. [17]. The unequal time correlator of the source 
is discussed in the next section. 

3. Generation of gravitational waves 

Gravitational waves are sourced by fluctuations in the energy momentum tensor which have a 
non-vanishing spin-2 contribution. The perturbed metric is of the form 

ds 2 = a 2 (-dr] 2 + (7^ + 2h ij )dx i dx j ) , (9) 

where 7^ is the metric of a 3-space of constant curvature K and hij is transverse and traceless. 
In Fourier space k l hij = h\ = 0. 

Einstein's equation to first order in gives (we set K = in what follows) 

(d 2 + 2Hd v + k 2 ) = MGa 2 ^ . (10) 

Here ITjj(k, 77) is the Fourier component of the tensor type (spin-2) anisotropic stress of the 
source and % = ^. A prime denotes the derivative w.r.t. conformal time 77. 

During a first order phase transition, anisotropic stresses can be generated by colliding bubbles 
and by inhomogeneities in the cosmic fluid (e.g. turbulence, magnetic fields) [16]. 

For a duration of the phase transition given by A77*, the final bubble size, i.e. the typical size 
of the bubbles when they start to coalesce and the phase transition terminates, is i?* = i^A^*, 
where Vb is the bubble velocity. Typically A77* < Hr x . As we shall see, depending on the 
unequal time correlation properties of the anisotropic stress, the peak frequency can be given 
either by the time scale, fc* ~ (A?7*) _1 , or by the spatial scale k* ~ i?" 1 . For the most promising 
case of a strongly first order phase transition, the bubble velocity is however close to the speed 
of light, Vb ~ 1 so that this distinction is not very relevant. 

4. The spectrum 

Because of causality, the correlator (11^- (771 , x)]n^ m (r/2, y)) = M.iji m (r)i,rf2,x — y) is a function 
of compact support. For distances |x — y| > max(?7i, 772), M. = 0. Therefore, the spatial Fourier 
transform, Miji m (r)i, 772, k) is analytic in k around k = 0. We decompose 11^ into two helicity 
modes which we assume to be uncorrelated so that the spectrum is even under parity 

Ilijfak) = egn+fak) + %n_(r?,k) 

(n + (r ? ,k)n;(r ? ',k')) = <n_(^,k)ni(»/,k')) = (27rfd 3 (k-k') P 2 x p Il (v,v / ,k) 

(n + (r ? ,k)n*_(r ? / ,k')) = o. 

Here px is the energy density of the component X with anisotropic stress n which has been 
factorized in order to keep k 3 Pn{i],rf ,k) dimensionless. Causality implies that the function 
Pu(r),r)' , k) is analytic in k. We therefore expect it to start out as white noise and to decay 
beyond a certain correlation scale fc*(?7, 77') > min(l/r/, 1/77'). 

If the gravitational wave source is active only for a short duration A77* (less than one Hubble 
time), we can neglect the damping term 2% in the equation of motion for h during the time 
when n is active. Here h denotes either of the two helicities of the gravitational wave and n is 
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Figure 1. The equal time anisotropic stress power spectrum from turbulence. On large scales 
the spectrum is scale invariant, on scales much smaller than the correlation scale it falls of 
like /c -11 / 3 which is determined by the Kolmogorov turbulence assumed for the velocity power 
spectrum. The black line is the numerical result and the red dashed line a simple fit. The 
normalization is arbitrary. Figure from [17]. 



the anisotropic stress source of the corresponding helicities. The solution with vanishing initial 
conditions is then 
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The gravitational wave energy density is given by 
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If the Universe is radiation dominated at the time 77* when the gravitational waves are generated, 
this leads to the gravitational wave energy spectrum 
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Here 



P n (uj,uj',k)= drj dv'PnMiky^-u'v'), (12) 

J— 00 J— 00 

On large scales, k < k* the GW energy density from a 'causal' source always scales like k 3 . This 
remains valid also for long duration sources. The only change for long lived sources is that the 
Green function e lk ^ v ~ v ' of the wave operator in Minkowski spacetime has to be replaced by the 
Green function of the wave operator in an expanding universe which in the case of a radiation 
dominated universe is given in terms of spherical Bessel functions. 

The scale l//c* is the correlation scale which is smaller than the co- moving Hubble scale 
l/k* < l/H* = The behavior or the spectrum close to the peak and its decay rate on 



smaller scales depend on the source characteristics, on its temporal behavior and on its power 
spectrum. 

Let us now briefly discuss the different results which have been obtained for the position of 
the peak of the gravitational wave spectrum. As we see from eqs. (12) and (11), to determine 
the the gravitational wave energy density we must know the unequal time correlator of the of 
the anisotropic stresses, Pn(i], rj' , k). This is not well known from numerical simulations but 
there exist different simple possibilities: 

• A totally coherent source is one for which the time dependence is deterministic. In this case 
the randomness is solely in the initial conditions and n(r/, k) = f(n, /c)n(r/„ I , k) for some 
deterministic function /. For the unequal time correlator Pr this leads to 

It has been shown in Refs. [18, 17] that the peak position of the GW spectrum oc 
k 3 Pn(k,k,k) for a totally coherent source is given by the peak of the temporal Fourier 
transform of the source, i.e. fc* = (A??*) -1 . 

• For a source with finite coherence time we can model the unequal time correlator of the of 
the anisotropic stresses, by 

Pn{rj,rj',k) = \/Pn(rf,rf,k)P n {'n,r],k)Q(x c - |r? - rf\k) , x c ~ 1 . 

This ansatz assumes that the source is totally coherent if the time difference is smaller than 
about one wavelength and uncorrelated for larger time differences. This ansatz resembles the 
suggestion by Kraichnan [19] for turbulence and has therefore been used for both turbulence 
and magnetic fields in Ref. [17]. However, this ansatz, like also the Kraichan decorrelation 
ansatz still lacks numerical confirmation. 

In this case, the peak of the GW spectrum is determined by the peak of the spatial Fourier 
transform of the source. 

• A totally incoherent source is uncorrelated for different times. Its unequal time correlator 
can be modeled as 

Pn(v, v', k ) = Pn(v', v', - rf) . 

Here the factor At/* is present to make the dimensions right, it can be replaced by some 
other short correlation time. In this case, the GW spectrum is also determined by the peak 
of the spatial Fourier transform of the source. 

The resulting gravitational wave spectra from the different hypotheses for the unequal time 
correlator are shown in Fig. 2. They are derived in detail in Ref. [17]. 

5. The electroweak phase transition 

According to the standard model, the electroweak transition is not even second order, but only 
a cross-over. Then, this transition does not lead to the formation of gravitational waves. Even 
though the Higgs vacuum expectation value does cross-over from zero to a finite value, the 
corresponding Higgs field fluctuations have a very short correlation length of the order of the 
inverse temperature and do not invoke significant gravitational field fluctuations. 

However, if the standard model is somewhat modified e.g. in the Higgs sector, or in certain 
regions of the MSSM parameter space, the electroweak phase transition can become first order, 
even strongly first order [7]. Then it proceeds via the formation of large bubbles of true vacuum. 
The bubbles themselves are spherically symmetric and do not generate gravitational waves but 
the collision events do. 




Figure 2. The GW energy density spectrum for a totally incoherent source (red, long-dashed), 
a totally coherent source (blue solid) and for a source with finite coherence time (black, short- 
dashed). The parameters are: % = 100 GeV, Arj*H* = 0.01, SlxMad = 2/9 ((v 2 ) = 1/3). The 
horizontal variable is K* = kR*/(2ir). Note that the peak position is roughly at k = ir/Ar]* = 
a few, hence A* ~ i?*/Ar/* = in the totally coherent case and at k = n/R* = a few, hence 
A'* ~ 1 in the incoherent case and in the case with finite coherence length. In order to clearly 
separate the peak positions, Vb — 0.1 has been chosen. Figure from [17]. 

In addition, the Reynolds number of the cosmic plasma is very high and the bubbles are 
expected to induce turbulence and, if there are small magnetic seed fields these are amplified 
by turbulence leading to an MHD turbulent plasma with turbulent kinetic energy and magnetic 
fields in equipartition [20]. 

The correlation scale is expected to be [21] of the order of ft* = j3 = (At?*) -1 ~ IOO/77* ~ 
10 _3 Hz, which is close to the frequency of the peak sensitivity for the space born gravitational 
wave antenna LISA, proposed for launch in 2018 by ESA [8]. 

5.1. Bubble collisions 

New numerical simulations of gravitational waves from bubble collisions during the electroweak 
phase transition [22] have given a slower decay law on small scales, like oc ft , see Fig. 3. 

In Ref. [18] it has been shown that the decays law ^f^r%r oc A; -1 is generically expected on 
small scales (large ft) if Pn{n,r],k) is continuous in time but not differentiable. This behavior 
seems reasonable for bubble collisions which start from vanishing overlap continuously but with 
a finite slope. 

5.2. Turbulent MHD 

Since the Reynolds number Re, of the cosmic plasma at T ~ lOOGeV is very high, the rapidly 
expanding bubbles of the broken phase provoke the formation of turbulence. Furthermore, in the 
broken phase the electromagnetic field does generically not vanish. Since the magnetic Reynolds 
number R m is even higher than the Reynolds number of the plasma [17], and thus the Prandl 
number [20] P m = R m /Re S> 1, the high conductivity rapidly damps the electric fields and we 
are left with a magnetic field in a turbulent plasma, MHD turbulence [20]. 
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Figure 3. d ^\^P from colliding bubbles, / = k/2ir, numerical results for f^x/^rad = 0.03. On 
the left panel the old results from Ref. [16] are shown while the right panel shows the result of 
Ref. [22]. Figure from [22]. 



To determine the anisotropic stress tensors from the turbulent flow and from the magnetic 
field which is the source of gravitational waves, we first discuss the vorticity and magnetic field 
spectra. Because both, the vorticity and the magnetic field are divergence free, causality requires 
that their spectra behave like k 2 for small k: 

(^(k)^(k')) = (27r) 3 5 3 (k-k')(k J k t -5 lJ )P v (k), (13) 
(BiMBjik')) = (27r) 3 5 3 (k-k')(k J k l -5 t ,)P B (k). (14) 

The functions (kjki — 5ij)P,(k) must be analytic at k = because of causality (the correlation 
functions are functions with compact support, they vanish for separations larger than the Hubble 
scales, thus their Fourier transforms must be analytic) hence P v (k) and Pe(fc) oc k 2 for small k, 
see [23] . Turbulence is steered during the phase transition and then decays freely with decay 
time of the order of A??*. This source is not truly long-lived, but since it decays only like a 
power law and not exponentially, it does still contribute significantly for r] > rj* + Ar/* and we 
expect a suppression factor e in eq. (8) which is larger than the naive factor (H*Ar]*) 2 . 

The behavior of the spectrum on scales smaller than the correlation scale k > fc* is expected 
to be a Kolmogorov spectrum for the vorticity field, P v oc k~ 11 / 3 and an Iroshnikov-Kraichnan 
spectrum [24] for the magnetic field, Pb oc k^ 1 1 2 . For the induced GW spectrum this yields 

dn GW .(k,f] ) ( njjj*) \ 2 j (k/h) 3 for k < k* 

d\n(k) ~ rad[r]0) \n r M J X \(k/h)- a for 

For • = v we have a = 11/3 — 1 = 8/3 and for • = B we have a = 7/2 — 1 = 5/2, see [25]. To find 
a we use that on small scale we expect the anisotropic stress spectrum to decay like the velocity, 
respectively magnetic field spectrum. With (|A| 2 ) oc Pu/k 2 and oc k 3 (\h\ 2 ) oc kPn so 

that a = 7 — 1, where 7 denotes the decay exponent of the anisotropic stress spectrum. 

In Figs. 4 and 5 the gravitational spectrum obtained from the wlectroweak phase trnasityion 
and from a first order phase transition at T* = 5 x 10 6 GeV are shown and compared with future 
gravitational wave experiments like LISA [8], the Bib Bang Observer (BBO) [26] asnd advances 
LIGO [14]. 
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Figure 4. The gravitational wave energy density from magnetic fields (red) and turbulence 
(blue, the sum is shown in black) generated during the electroweak phase transition at 
T* = lOOGeV. The time-decorrelation of the source is modeled by the finite correlation approach 
(Kraichnan decorrelation) . The source energy density spectrum is assumed to be maximal, 
^W^rad = 2/9, see [17]. Equipartition in turbulent kinetic energy and magnetic field energy is 
assumed. The sensitivity curves are from [27]. Figure taken from [17]. 



5.3. Limits of primordial magnetic fields from the EW transition 

It is difficult to estimate Qb(v*) or ^v(v*) directly, but since causality requires the spectra to be 
so blue, the limit on gravitational waves which comes from small scales k — k*, yields very strong 
limits on primordial magnetic fields on large scales already from the simple nucleosynthesis 
constraint, Q gw <, 0. lSl ra( j . If the magnetic field spectrum behaves like Pb(/s) oc k n , the 

gravitational wave energy density per log interval must go like ^^^j* - oc fc 3+n on large scales. 
For causal magnetic fields we have n = 2, hence the magnetic field energy spectrum is very blue, 

^l^ cxfc 5 for k<K. (15) 

The main contribution to the gravitational wave energy density comes from the correlation 
scale, Q gw ~ d d\n(k) " > wri ereas we are interested in the magnetic fields on much larger scales, 
for example = (O.lMpc) -1 . Using our simple formula (8) with e ~ 1 the nucleosynthesis 

limit yields k^ 2 B(k*, rjo) < 10~ 6 Gauss. The time 770 denotes today. For the electroweak phase 
transition A;* ~ 100%* — 10 _3 sec _1 while we are interested in magnetic fields on scales of O.lMpc 
giving ki = 10(Mpc) _1 ~ 10 _13 sec _1 . With the above scaling, the limit on a magnetic field on 
this scale is 

( k \ 5//2 

k 3 L /2 B{k L ,ri ) < ( -± ) 10" 6 Gauss = 10" 31 Gauss . (16) 



Note also that this limit is very stable with respect to the small unknown pre-factor e. Since 
VL gw oc .B 4 , changing e by 2 orders of magnitude yields a change of only a factor of 3 in the limit 




Figure 5. To the GW background from the electroweak phase transition also shown in Fig. 4 
(here left black line) we add the result for a phase transition at T = 5 x 10 6 GeV with 
A 77*%* = 0.02 (black curve to the right). In this case, the gravitational wave background 
might even lie within the reach of advanced LIGO. The sensitivity curves are from [27]. 



for B. 

This limit is unavoidable if the magnetic field on large scales changes simply by flux 
conservation so that B(rj) = B(rj*)(a*/a) 2 . Scaling this to today yields B(r]o) = -B(r/*)/(l + £*) 2 - 

The limit is so severe because the magnetic field spectrum is so red. The simple requirement 
that pB(t*) < Prad(^*) yields roughly the same limit. However, this limit is not entirely save 
since at very early time most of the energy density might in principle have been in the form of 
small scale magnetic fields which later are converted into heat by viscosity damping on small 
scales. The fact that such magnetic fields would generate a sizable background of gravitational 
waves, however, forbids this possibility. This has been realized for the first time in Ref. [28]. 
The inferred limits on magnetic fields as a function of the spectral index are shown in Fig. 6. 

5.4- Helicity 

The most optimistic dynamo models for the generation of the magnetic fields of about 10 -6 Gauss 
which are observed in galaxies and clusters, need seed fields of the order of at least 10 -21 Gauss. 
These are excluded by the above limit. There is a possible way out: If the magnetic field has 
non-vanishing helicity it can develop an inverse cascade and power can be transferred from small 
to larger scales [29]. In this case the large scale magnetic fields do not simply decay by flux 
conservation in an expanding universe but they can feed from the small scale fields by inverse 
cascade. 

During the electroweak phase transition parity is broken. Actually, the Chern-Simons 
winding number of the gauge field, Ncs °c / F A A, which is related to the baryon number, 
has an electromagnetic part to it which is simply the helicity of the magnetic field [30], 
H = V^ 1 f v (A ■ B)d 3 x . It is interesting to note that this relates the baryon number and 
the magnetic helicity. It has been shown that the parity violation in the stress energy tensor 
coming from such helical magnetic fields lead to T-B and E-B correlations in the temperature 
and polarization spectra of the cosmic microwave background [31]. By the same token they also 
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Figure 6. The nucleosynthesis limit on magnetic fields from pB{t nuc ) < 0.1/9 ra d(i n uc) (black solid 
line), and from the gravitational wave constraint (red lines) as a function of the spectral index 
n. The limits for magnetic fields from the electroweak phase transition are short-dashed (red), 
while the long dashed (red) line corresponds to magnetic fields from inflation at a GUT scale. 
Causality requires that the spectral index be n = 2 for magnetic fields from the electroweak 
phase transition while n is in principle undetermined for inflation. Figure from [28]. 

generate gravitational waves with non- vanishing helicity [31]. 

The inverse cascade which is caused by helicity conservation for a helical field and is absent 
for non- helical magnetic fields has been investigated numerically, e.g., in [32] and is shown in 
Fig. 7 

In Ref. [32] simple power laws are derived for the evolution of the correlation scale in the 
inverse cascade regime. Using this simple prescription, the mitigation of the magnetic field limits 
by an inverse cascade of maximally helical magnetic fields has been estimated in Ref. [33]. It 
has been found, however, that the inverse cascade is not sufficient to present a way out for the 
magnetic fields from the electroweak phase transition, but it can work for helical magnetic fields 
generated during the QCD phase transition. The QCD phase transition might be first order if 
the neutrino chemical potential is large, as discussed in [10]. 

If the magnetic field is helical, the GW background would not be parity symmetric. There 
would be more GW's of one helicity than of the other [31]. 

6. Conclusions 

• The rapid expansion and collision of bubbles during a first order phase transition stirs the 
relativistic cosmic plasma sufficiently to lead to the generation of a (possibly observable) 
stochastic gravitational wave background. 

• Observing such a background would open a new window to the early Universe and to high 
energy physics. 
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Figure 7. Left: the MHD evolution of the large scale magnetic field energy density spectrum 
with vanishing helicity. The field is simply damped by viscosity on small scales. Large scales 
are not effected. The viscosity scale is growing with time. 

Right: the case of non- vanishing helicity. First, the helicity is relatively small and the evolution 
of the helical field is not very different from the non-helical case. But due to helicity conservation, 
helicity rapidly becomes maximal and the evolution enters the inverse cascade regime. 
For simplicity, the correlation scale is set equal to the damping scale. Figure from [32]. 



Generically, the density parameter of the GW background is of the order of 

tt gw (t ) ~ e $7 ra d(*o) . n 

where e < 1 is determined by the duration of the source, e ~ 1 if the decay time of the 
source is roughly a Hubble time. 

The spectrum grows like ^^(k)^ ^ ^ on l ar S e scales and decays on scales smaller than 
the correlations scale ~ I/A77*. The decay law depends on the physics of the source. 
Within the standard model, the electroweak phase transition is not of first order and does 
(probably) not generate an appreciable gravitational wave background. However, simple 
deviations from the standard model can make it first order. 

In this case we expect a GW background which can marginally be detected by the LISA 
satellite. 

It has been proposed that the magnetic fields generated during the electroweak phase 
transition, could represent the seeds for the fields observed in galaxies and clusters. 
Unfortunately this idea does not withstand detailed scrutiny. The magnetic fields cannot 
have enough power on large scale to represent such seed fields. On small scales the magnetic 
field power which is initially present is dissipated later on by the viscosity of the cosmic 
plasma. 

If there is no inverse cascade acting on the magnetic field spectrum, the limits on the large 
scale fields coming from the generated GW background are too strong to allow significant 
magnetic fields even for the most optimistic dynamo mechanism. 

If the magnetic field is helical, helicity conservation provokes an inverse cascade which 
can alleviate these limits. Most probably the relaxed limits are still too stringent for 
magnetic fields from the electroweak phase transition, but magnetic fields from the QCD 
phase transition might be sufficient to seed the fields in galaxies and clusters if they are 
helical. 

Helical magnetic fields generate a parity violating gravitational wave background, 
\h + (k)\ 2 / |/i_(/c)| 2 . This parity violation is in principle observable. 
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